Tailorable couplings of a cantilever with a superconducting charge qubit: 

Quantum state engineering 
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We propose a theoretical scheme to realize tailorable couplings between a cantilever and a superconducting 
charge qubit. By tuning the controllable parameters of the qubit, both linear and nonlinear couplings between 
the cantilever and the qubit can be achieved. Based on these couplings, we show the preparation of the cantilever 
into some interesting quantum states, such as superposed coherent states and squeezed states, via manipulating 
and detecting the qubit. We also study the influence of the environment on quantum states of the cantilever. 
It is indicated that decoherence induced by the environment can drive the cantilever from superposed coherent 
states into the steady coherent state. It is also found that the environment can induce the steady-state position 
squeezing of the cantilever under a critical temperature. These results will shed new light on production of 
nonclassical effects of the cantilever. 

PACS numbers: 85.85. +j, 85.25.Cp, 42.50.Dv 



I. INTRODUCTION 



As is well known, it is of very significance to couple a me- 
chanical object to an electronic system since the electronics 
may be used to measure the quantum nature of the mechan- 
ical object while quantum effects in electronic systems may 
be measured by mechanics. Cantilevers, as important compo- 
nents of magnetic resonance force microscopy (MRFM) fljj- 
3], have attracted much attention of both theorists and experi- 
mentalists in many fields of physics such as condensed matter 
physics and quantum information science |4j-|6j]. For exam- 
ple, MRFM has been recently proposed as a qubit readout de- 
vice for spin-based quantum computers 10, [H]. As one kind of 
nanomechanical resonators, the cantilever can also be used as 
a platform to study some fundamental problems in quantum 
mechanics such as quantum measurement jfj, quantum de- 
coherence [ 10], the boundary between quantum and classical 
communities lITll] . and test of quantum mechanics in macro- 
scopic scales lll2l [l3ll . By far, with the high-speed develop- 
ment of modern micro-fabrication techniques, the preparation 
of nanomechanical resonators with high frequency and high 
quality factor lfl4l - [l6ll has become pos sible. At the same 
time, efficient cooling technologies ltT7l430ll drive nanome- 
chanical resonators approaching quantum realm. Therefore, 
it is expected to observe quantum evidences [31] in nanome- 
chanical resonators such as cantilevers. As a precondition, 
it is an interesting topic that how to generate some nonclas- 
sical states Il32l - l45ll such as Fock states, superposed states, 
squeezed states, and entangled states in cantilevers. In addi- 
tion, from the viewpoint of MRFM, cantilevers play the role 
of probers to read out the states of a single spin, therefore 
the preparation of cantilevers into some states with position 
squeezing can improve the measurement precision 14611 . 

For manipulation of the vibrational states of nanomechan- 
ical resonators, many schemes have been proposed by cou- 
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pling nanomechanical resonators with vario us phys ical sys- 
tems such as superconducting qubits IU2U471450I1 . super- 
conducting transmission line resonators Il5lij53ll . cold ions, 
atoms and molecules ll54l - l57tl . and quantum dots ll33l l58l - 
l60ll . These considered physical systems play the role of in- 
directly assistant controller 16111 . Therefore, they should be 
of well controllability and readability. For example, during 
the latest decade, great advances in quantum information pro- 
cessingbased on superconducting charge qubits have been 
made 1 62, 63]. It has been shown that superconducting charge 
qubits have well controllability and readability lf64Tl . Stimu- 
lated by these, in this paper, we propose a new scheme to cou- 
ple a cantilever with a superconducting charge qubit lr65ll . We 
can manipulate the cantilever through tuning and measuring 
the qubit as an indirect controller. Both linear and nonlinear 
couplings ll39ll between the cantilever and the qubit can be ob- 
tained. Especially, we emphasize that the obtained nonlinear- 
type coupling for a cantilever is one of the main results in this 
work. Based on these couplings, we show how to create su- 
perposed coherent states and squeezed states of the cantilever 
by manipulation of the qubit. The preparation of superposed 
coherent states shows the appearance of quantum superposi- 
tion in the cantilever. And the squeezing in the cantilever not 
only shows nonclassical evidence but also has wide potential 
for practical application. 



This paper is organized as follows. In Sec. [II] we intro- 
duce the physical model and present correspondent Hamilto- 
nian. In Sec. [TTTJ we show how to generate superposed co- 
herent states of the cantilever, and investigate the influence of 
decoherence induced by environment on superposed coherent 
states of the cantilever. In Sec. IIV1 we show how to obtain 
a dynamical squeezing of the cantilever, especially it is indi- 
cated that there exists a steady-state position squeezing. We 
conclude this paper with some discussions in the last section. 
Finally, we present two Appendixes|A]and|B]for derivation of 
evolution equations of the cantilever in an environment. 
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II. PHYSICAL MODEL AND HAMILTONIAN 

We start with introducing the physical setup as illustrated 
in Fig. [T| a cantilever is fabricated above a superconducting 
charge qubit which is formed by an SQUID-based Cooper- 
pair box l65ll . A ferromagnetic particle mounted on the can- 



NAMR * ! 

X-/ 

* i Magnetic tip 




FIG. 1 : (Color online) Setup of a cantilever mounted with a magnetic 
tip coupling with an SQUID-based charge qubit. The magnetic field 
generated by the magnetic tip threads through the superconducting 
loop of the charge qubit. Since the generated magnetic field depends 
on the vibration amplitude z of the cantilever, and the magnetic flux 
in the loop is a controllable parameter of the qubit, therefore a cou- 
pling between the cantilever and the qubit is induced. 

tilever tip produces a magnetic field i6r3l 

fitip = ^3 [3 (ft • m) ft-rn], (1) 

threading the superconducting loop of the qubit, where no is 
the vacuum magnetic conductance, « is the unit vector point- 
ing to the direction from the tip to the center of the loop. Here 
we assume that the magnetic tip is right on top of the cen- 
ter of the loop. The vector in is the magnetic moment of the 
magnetic tip pointing to the z direction, and r is the distance 
between the tip and the center of the loop. For a tiny vibration 
of the cantilever, the magnetic field threading the loop of the 
qubit can be approximated as 

fi tip » B - Cz, (2) 

where Bq = fj.o'n/(2nr 3 ) and C - 3fiom/(2nr 4 ) with m = \m\. 
Since the generated magnetic field depends on the position z, 
then couplings between the cantilever and the qubit can be 
induced as follows: The vibration of the cantilever leads to a 
change of the magnitude for the magnetic field threading the 
loop of the qubit. As a result, the corresponding magnetic flux 
will be changed. At the same time, since the magnetic flux is 
a controllable parameter of the qubit, therefore the vibration 
can induce a coupling between the cantilever and the qubit. 

As for the SQUID-based charge qubit [65], its Hamiltonian 
reads 

H g = \E C \n g - cr z - Ej cos °"*> (3) 

where Ec = e 2 /[2(2Cy + C g )] is the single-electron charg- 
ing energy, Cj and C g are, respectively, the capacitances of 



each Josephson junction and the gate capacitance in the qubit; 
n g = C g V g /2e is the gate charge number with V g being gate 
voltage; Ej is the Josephson coupling energy; <f> e and 3>n are 
the externally biasing flux and flux quanta, respectively. The 
Pauli operators introduced in Eq. ([3]) are defined by 

<r z = 10X01 — 11><1|, ^, = 10X11 + 11X01, (4) 

where the states |0) and |1) represent that there is no and one 
extra cooper pair on the inland, respectively. 

We consider a case that the externally biasing magnetic flux 
O e is composed of two parts. One is generated by the mag- 
netic tip, and the other is generated by an externally control- 
lable electric current. We can express the total biasing mag- 
netic flux as (D e = (Z?ti p + B X )S , where B x is the magnetic field 
generated by the externally controllable electric current, S is 
the area of the superconduting loop. For a tiny vibration, the 
cantilever can be modeled as a quantum harmonic oscillator, 
which is depicted by the usual Bosonic creation and annihi- 
lation operators a ' and a, satisfying the commutative relation 
[a, a'] = 1. Then the magnetic field generated by the magnetic 
tip can be expressed as 

Btip ~ B - Czo(a T + a), (5) 

where zo = 1/ ^2pm) (with fi = 1) is zero-point uncertainty 
for the ground state of the cantilever, fj. and a> are the mass 
and frequency of the cantilever, respectively. 

The Hamiltonian of the total system including the cantilever 
and the charge qubit reads 

H — ~^~ cr z + a>a ' fl ~ Ej cos[0o + cf)(a ' + a)]cr x , (6) 

where a>o = 8Ec(n g - 1/2) and co are, respectively, the fre- 
quencies of the qubit and the cantilever. We also introduce 
two parameters 

0o = ttS (Bq + JJjO/Go, = -nS Czo/Oo- (7) 

Hamiltonian © obviously shows a nonlinear coupling be- 
tween the cantilever and the charge qubit. A coupling of 
similar form between a bosonic mode and a two-level sys- 
tem has been obtained in a trapped-ion system i67tl . A re- 
cent scheme has been proposed to obtain a nonlinear interac- 
tion between a doubly-clamped beam and a superconducting 
charge qubit 113911 . However, the method proposed in Ref. fl39ll 
is not valid for a cantilever. 

Hamiltonian <JSJ is very useful in quantum information pro- 
cessing. Many useful interactions can be tailored from Eq. © 
by choosing proper parameters. For example, we tune the ex- 
ternally controllable current such that cos <po - 0, then Eq. (O 
becomes 

t ^o t 
Hi — Loa'a + ~2°~z + Ejsin[cf)(a' + a)]cr x , (8) 

which can be reduced to the well-known Jaynes-Cummings 
Hamiltonian without rotating wave approximation 

Hi » -^-cr z + dLta'fl + g(a' + a)cr x , (9) 
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by expanding the sine function up to the first order of param- 
eter <p, where we introduce the coupling strength g — Ej<p. 

On the other hand, if we choose the external magnetic flux 
to ensure cos 0o = T then Eq. © reduces to 

H2 — toa'a + -T-°"z - -Eycos[0(a' + a)]cr x , (10) 

which can be further simplified by expanding the cosine func- 
tion up to the second order of <p, 

H2 ~ cja'a + ~^°~z ~ Ejcr x — g (a' + a) cr x , (11) 

where g' = Ej(/> 2 /2 is a nonlinear coupling strength between 
the cantilever and the qubit. These two kinds of couplings 
given in Eqs. (0 and (ITTb are very useful in quantum optics 
and quantum information processing. As examples, in the fol- 
lowing two sections, we will study quantum state engineering 
based on these couplings. 



III. PREPARATION OF SUPERPOSED COHERENT 
STATES 

Superposed coherent states are typical quantum states 
which exhibit nonclassical properties. In this section, we show 
how to prepare the cantilever into superposed coherent states 
with the above obtained Hamiltonian We also study the 
decoherence of the created superposed coherent states when 
the cantilever is subjected to an environment. 



A. Generation of superposed coherent state 

Firstly, we consider an ideal situation in which there is no 
dissipation for the cantilever. Since the state preparation can 
be realized in a very short time interval, it is reasonable to 
neglect the dissipation during the state preparation process. 
We tune the gate voltage V g such that n g = 1/2, that is a>o = 0, 
then Eq. (0 reduces to the conditional displacement harmonic 
oscillator (CDHO) Hamiltonian 



where we introduce the parameters 



^cdho - a>a'a + g(a' + a)cr x . 



(12) 



Corresponding to the qubit in states |+), the displacement 
terms are ±g(a' + a), respectively, where states |±) are the 
eigenstates of Pauli operator <r x , with the respective eigenval- 
ues ±1. 

For generation of superposed coherent states, we suppose 
the total system consisting of the cantilever and the qubit is 
initially prepared in a state \<p(0)) = |a,) ® |0) ? , where |a,) 
is the usual Glauber coherent state, which is defined as the 
eigenstate of annihilation operator a, i.e., a\ai) = a,[a,), and 
|0> is defined by |0> = (|+> + |-»/ V2. Making use of Hamil- 
tonian (fT2l . the state of the total system at time t is 

\<p(f)) = exp(-;r#cDHo)W0)> 



a + ) +e w -\a. 



\e w *\a + )- 



>)|0> 
ar_>) |1>] , 



a ± = a ie -°" ± ^(e^' - 1), 

(!) 

0± = - \gt~(-±ai)sm(a>f) 
to I \oj / 



(14a) 
(14b) 



During the derivation of Eq. ( fT3l . we have used the follow- 
ing formula M68I1 . 



,[6»03ia+/3 2 a t a+Aa t )] 



g/ifl 1 e fia r a e f 3 a e fi 



with 









h 






h 




-i)/a. 


n 







(15) 

(16a) 
(16b) 
(16c) 
(16d) 

The action of the operator given in Eq. ( fT5l l on a coherent state 
la,) yields jH: 

e im^a + fiW)] la . } = e , + a . e h} t (17) 
where s is given by the following expression 

e = / 4 +/3tf/ + (|a;exp(/ 2 )| 2 -N 2 + |/i| 2 

+2Re[/ia;exp(/ 2 *)])/2. (18) 

From Eq. ( fT3l , it is obvious to prepare the cantilever into 
superposed coherent states through measuring the qubit. Cor- 
responding to the states |0) and |1) of the qubit are measured, 
the cantilever collapses to the following superposed coherent 
states 



\ilf ± ) = N ± {e ie +\a + )±e ie -\tt-)).. 



(19) 



where N~ 2 = 2{1 + Re[exp(-/(6» + - 0_))<ar+|a_>]} are normal- 
ization constants. 

For a special case, we suppose the cantilever is initially pre- 
pared in a vacuum state, i.e., a, = 0, then the states of the can- 
tilever at time t are the so-called Schrodinger cat states l69ll 



\<P±)=M ± (\/3): 



(20) 



(13) 



withyg = g[exp(-/«f) - l]/co and Mr 2 = 2[1 + exp(-2jS 2 )]. 

When the qubit is detected in states |0) or |1), the coupling 
term g(a f + a)cr x between the qubit and the cantilever will en- 
tangle them. Therefore, from the experimental viewpoint, we 
should decouple the cantilever and the qubit, as long as the 
superposed coherent states are prepared. The method to de- 
coupling is making the second measurement on the qubit in 
states |±). Since states |±) are eigenstates of the operator <x v , 
then the qubit will stay in states |+) forever, and the dynamics 
of the cantilever is governed by the displaced harmonic oscil- 
lator (DHO) Hamiltonian 

#dho = wa'fl +g s {a ] +fl), (21) 

where g s - +g corresponding to the qubit in states [+). 
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B. Decoherence of superposed coherent states 

In this subsection, we investigate the decoherence of super- 
posed coherent states produced in the previous subsection. As 
a practical physical system, the cantilever couples inevitably 
with its external environment. Therefore, the cantilever pre- 
pared in superposed coherent states will loss its coherence and 
energy. We suppose that the preparation time for this initial 
state is very short, thus we neglect the decoherence in the 
course of the initial state preparation process. 

The dynamics of the cantilever is governed by the quantum 
master equation 

P = i\p,H DHO ]+£[p], (22) 

where the decoherence of the cantilever is phenomenologi- 
cally represented by the superoperator £. At a temperature of 
T, this superoperator can be written as 04611 

JL\p\ — —{n,h + \){2apa—a'ap—pa'a) 

J - t t 
+— n / /,(2fl'pa — aa' p — paa' ), (23) 

where y is decay rate, and n,/, = l/[exp(a>/T) - 1] is average 
thermal excitation number of the thermal bath at frequency at. 
Equation d22b shows that there are three kinds of actions on 
the cantilever. The free Hamiltonian oxva rotates the system 
in phase space. The driving term g s (a ' + a) displaces the can- 
tilever in phase space. And the dissipation term X decreases 
the coherence and energy of the system. 

To see the evolution of the cantilever, we need to solve 
quantum master equation d22i i. In order to do this, we intro- 
duce the following transform, 

p {1 \t) = D\a{t))R\6(t))p{t)R{6{t))D{a(t)), (24a) 

R(6(i)) = exp(-!0(r)« + tf), (24b) 

D(a(t)) = exp(a(f)fl f - a*(t)a), (24c) 

with 

0{t) = cot, (25a) 

a(t) = -^-(e^ -e**), (25b) 

then quantum master equation d22b can be transformed to a 
standard form 

fP = Z(n + l)(2ap (2) fl t - atapW-pfWd) 

+ |n(2a t p (2) fl - aa'p {2) - p (2 W), (26) 

which describes the evolution for a harmonic oscillator in a 
heat bath. The detailed derivation of quantum master equa- 
tion ( f26b will be presented in AppendixlAl 

The operators R(6) and D(a) in Eq. d24l i are, respectively, 
the usual rotation and displacement operators for a harmonic 
oscillator in phase space. In principle, the solution for quan- 
tum master equation d26l > can be obtained with the superoper- 
ator method. However, for simplicity, we only give the analyt- 
ical solutions for the zero temperature case in the following. 



Based on the above discussions, the dynamical evolution of 
the cantilever can be obtained as follows: When the cantilever 
is initially prepared in an initial state p(0), then the state p(f) of 
the cantilever at time t can be obtained through the processes 

p(0)^p (2) (0)^p (2) (?)^p(0. (27) 

The two processes p(0) — > p^ 2) (0) and p <2) (?) — > p(f) are deter- 
mined by the transform given in Eq. ( 124b and its inverse trans- 
form, and the evolution process p (2) (0) — > p (2) (f) is governed 
by quantum master equation ( f26b in the transformed represen- 
tation. 

We assume that the initial state of the cantilever is 

|^(0)> = %(^> + e ''n-/3>), (28) 

where N v = [2(1 + exp(-2y6 2 ) cos {f))Y 1 ^ 2 is the normaliza- 
tion constant. Then the state of the cantilever at time t is 

p{t) = M 2 (\/3 + ){/3 + \ + |/L></?-| + exp(A r + x^W+WA 

+^exp(A r -jA ! )^->^ + |), (29) 

where the parameters are defined as 

A r = 2\fl\\e-y< - 1), (30a) 
Ai = 2lm{a{t)l3*e-^), (30b) 
P ± = (a(t) ±pV?)e _i< <*. (30c) 

The detailed derivation from state d28l > to state d29l ) will 
be given in Appendix [B] From Eq. d29l ), we can see that 
p\(oo) = y3_(oo) = —igs/(y/2 + icj) at the long time limit. 
Therefore the steady state of the cantilever is a coherent state 
\-ig s /(y/2 + iai)), which is resulted from the net actions of 
the free evolution, the coherent driving, and the decoherence. 
When g s = 0, the steady state reduces to vacuum state |0), 
which implies the cantilever approaching an equilibrium with 
the zero temperature environment. 

For seeing clearly the decoherence of the generated super- 
posed coherent states, we address the evolution of the Wigner 
function for these states. The definition of the Wigner func- 
tion H70n of a density operator p is 

W(g) = 2Tr [D(-£)pD(£)e' Wfl ] . (31) 

For the state given in Eq. d29l . the Wigner function is obtained 

as 

W(g) = 2JV 2 (exp(-2|£-/? + | 2 ) + exp(-2|^- J S_| 2 ) 

+2Re [e- iip exp(A r + /A,- + ©)]) , (32) 

where is given by the following expression 

= 2 (#T + f p\ - |^| 2 ) - I (^6 + | 2 + \/3-\ 2 + 2/T p\) . (33) 

In Fig. [2] we plot the Wigner function given in Eq. (f32b at 
different times for the case of <p = and g s = g. When t = 0, 
the Wigner function is right for the state \tp + ) given in Eq. (l20b 




FIG. 2: (Color online) Plot of the Wigner function W(^) given in Eq. J32t for the state of the cantilever at different scaled times r = yt: (a) 
T = 0, (b) t = 0.06, (c) t = 0.5, (d) r = 0.1, (e) t = 1.5, (f) t = 20. Other parameters are set as co/y = 100, g/y = -300, /3 = 3, and ip = 0. 



with yS = 3. This state is the so-called even Schrodinger cat 
state. We can see from Fig. 0a) that there is some coherence 
in state ( 120) . This is because the Wigner function in Fig. [2J a ) 
exhibits some interference fringes (namely some oscillations), 
and these oscillations imply quantum coherence in the state. 
With the increase of the time t, we can find that the oscil- 
lations decrease gradually. At the same time, the positions 
of the two main peaks of the Wigner function rotate on the 
phase plain and move gradually to a point. This point corre- 
sponds approximately to the steady state of quantum master 
equation ( l22l i for T — 0. In fact, this steady state is a co- 
herent state | - igJ(y/2 + US)) ~ |2.999 + 0.015/), which is 
obtained from Eq. d29l by taking the long time limit. More- 
over, with the dissipative evolution, the negative values of the 
Wigner function will disappear gradually, which implies that 
the nonclassical properties of the cantilever decreases with the 
decoherence. Therefore, the actions of the environment and 
the driving force will destroy the coherence of the superposed 
coherent states and drive the cantilever into a steady coherent 



state. 



IV. DYNAMICAL SQUEEZING 

In the above section, we have study the creation of super- 
posed coherent states based on the obtained linear Hamilto- 
nian. In this section, we study the creation of dynamical 
squeezing as an application of the nonlinear Hamiltonian. The 
created squeezed state not only exhibits nonclassical proper- 
ties, but also is useful for precise measurement. 

A. Dynamical squeezing without dissipation 

From Hamiltonian (fTTT i. we control the gate voltage V g such 
that coo = 0, then Hamiltonian ( fTTT ) becomes 

H = ua'a - [Ej + g'(a< + a) 2 ] cr x . (34) 
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It can be seen from Hamiltonian ([34-b that if the qubit is ini- 
tially prepared in one of the two eigenstates |±) of o~ x , then 
the qubit will stay in this state forever. Corresponding to the 
cases of the qubit in |+), the conditional Hamiltonians of the 
cantilever are 



H k — a> k a^a + g k (a' 2 + a 2 ) , k = ±1, 



(35) 



where co k — to + 2g k and g k = -kx g'. 

To diagonalize Hamiltonian d35b . we introduce a unitary 
operator lT7ltl 

exp(-S k ) = exp [-A k (a 2 - a 1 " 2 )] . (36) 

Using the commutative relations [a, S k ] = -2A k a^ and 
[a^Sfc] = -2A k a, we obtain 

e~ Sk ae St = a cosh(2^) - a 1 " smh(2A k ), (37a) 
e~ Sk a f e Sk = a f cosh(2/t it ) - a sinh(2^). (37b) 

Application of the transform defined in Eq. (|36T > on Hamilto- 
nian (l35l) leads to 



(38) 



H k = exp(-S k )H k exp(S k ) 

= Q k a ] 'a + k k (a n + a 2 ) + C k , 

where we introduce the following parameters 

Cl k = co k cosh(4/U) - 2g k smh(4A k ), 
A k = -u> k sinh(4/} J t)/2 + g k cosh(4A k ), 
C k = u k sinh 2 (2ii) - g k sinh(4^). 



(39a) 
(39b) 
(39c) 



By choosing proper parameter A k to ensure = 0, namely 
tanh(4/ti) = 2g k /u> k , then we may choose 

sinh(4^) = cosh(4,l*) = ^, (40) 

Ll k Ll k 



where Q. k = ^u) 2 k - 4g 2 k . Therefore we obtain the diagonal- 
ized Hamiltonian 



H k = Q. k a'a, 



k = ±1. 



(41) 



The unitary evolution operator relating to Hamiltonian (|4"TT i is 
V k (t) = exp(-iH k t). Notice that in Eq. (PTTT i we have discarded 
the constant term C k . 

For investigation of the squeezing of the cantilever, we as- 
sume the cantilever is initially prepared in a coherent state \a) 
and the qubit in states |±), where the coherent amplitude a is 
assumed to be a real number for simplicity. After a coherent 
evolution of a time f, the cantilever evolves into state 



|<D t (f)> = e St V k {t)e- Sl \a). 



(42) 



Using the equation a \ a) - a\a) and \a) - e Sk VAi)e ^'lO^f)), 



we obtain 
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FIG. 3: (Color online) Plot the relative fluctuation {Az) 2 _ l lz^ given in 
Eq. J47I ) against the scaled evolution time ( = a>t. Other parameter 
are set as g'/a> = 0.0115. The relative fluctuation (&zj 1 _ l lz\ < 1 
means a squeezing of the position z for the cantilever. 



which can be further written as 

A k (f)\®k(f)) = a\® k (i)), (44) 

with 

A k (t) = e Sk V k (f)e- Sk ae Sl Vl(f)e- Sk 

= [cos(Q. k t) + i cosh(4A k ) sin(Q. k t)]a 

+i sinh(4/I t ) sin(Q /t f)a 1 ■ (45) 

The evolution of the cantilever follows the eigenstate of the 
quasi-excitation operator A k (t) with eigenvalue a. According 
to Eq. (1451 1. we can express the operators a and aJ in terms of 
A k (t) wd A\{t), 



a = [cos Q^f + / cosh(4/ljt) sin(Qjtf)]Ajt(f) 
+i sinh(4/U) sin(Qi;f)A](f). 



(46) 



Based on Eq. d46l ), the fluctuation (Az) 2 = (z 2 } k - {z) 2 k of the 
coordinator z for the cantilever is obtained as 



(Az) 2 (f) = z 2 [cosh(4A*) - sinh(4^)] 2 sin 2 (Q^f) 
+ cos 2 (QiO- 



(47) 



e Sk V k (t)e- Sk ae Sk vl(t)e- Sk \Q> k (t)) = ar|O*(0>, 



(43) 



From Eq. Wl\ , we find that, only in the case of A; = -1, 
the squeezing of the position z of the cantilever occurs. In 
Fig. [3] we plot the relative fluctuation (Az) 2 j/zq as a function 
of the scaled evolution time f = cot. Figure[3]shows a periodic 
squeezing with the scaled evolution time f . 



B. Dynamical squeezing with dissipation 

In the above subsection, we study the dynamical squeezing 
for the ideal case in which there is no dissipation. However, 
any systems will couple inevitably with the environment. In 
this subsection, we consider the dynamical squeezing of the 
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cantilever by taking the environment into account. In the pres- 
ence of an environment, the evolution of the cantilever is gov- 
erned by the following quantum master equation 



p = i[p(t),H k ]+£p(f), 



(48) 



where the Hamiltonian H k has been given in Eq. d35l l. and the 
superoperator has been given in Eq. d23l l. Based on quantum 
master equation (|48]>, we can obtain the following equations 
of motion, 



(a) = -i<o k (a) - 2ig k (a ] ) - -(a), 

t t y i 

{a 1 } = iu> k (a<} + 2ig k {a) - ^(a 1 ), 



dt 
d_ 

dt 
d 
dt 

^r{a n ) = -(y - 2iu) k ){a a ) + Aig^a) + 2ig k , 
dt 



(a ) = -{y + 2iu) k ){a > - Aig k (a { a) - 2ig k , 



(49a) 
(49b) 
(49c) 
(49d) 



— (a' a) = -yW'a) - 2ig k (a n ) + 2ig k (a 2 } + yh th . (49e) 
dt 

According to the definition of z - Zo(a + a ' ), the relative fluc- 
tuation of the coordinator operator z can be written as 



(MO), 2 



= [{a 2 ) - {a) 1 + (a n ) - (a 1 ") 2 
+2<fl 1 a> - 2<fl t )<a) + l) . 



(50) 



Therefore, under a given initial condition, we can obtain the 
solutions of Eq. d49l . and then the evolution of the fluctua- 
tion (Az) 2 (f) can be obtained. In the following we suppose the 
cantilever is initially prepared in a vacuum state |0). 

In principle, the solutions of Eq. d49l can be obtained. 
However, we do not present the solutions here since there are 
very complicate. In Fig. |H we plot the relative fluctuation 
(Az) 2 j(?)/Zq against the time t for different temperatures. It 
can be seen from Fig.|4]that the relative fluctuation (Az)\/zk 
evolves gradually approaching a steady-state value with the 
increase of the time t. For a short time, the relative fluctuation 
evolves with some oscillations. At the same time, the relative 
fluctuation increases with the increase of the bath temperature 
T. Therefore, we can obtain a conclusion that the high tem- 
perature T can destroy the squeezing for the position z of the 
cantilever. 

For obtaining the steady-state properties of the squeezing, 
we obtain the steady-state solution of Eq. d49l as 



<a(oo)> = (a'Xoo)) = 0, 

2g k {2n lh + l)(2a>i + iy) 



<« z (<*>)> 



<a |2 (oo)> = - 



r 2 +4w 2 -16g 2 ' 
2g k (2h, h + l)(2a) k - iy) 



y 2 +4a, 2 -16s 2 



(a'a(oo)) = n t h + 



8g 2 k (2ii th + 1) 
y 2 +4w 2 - 16g 2 



(51a) 
(51b) 

(51c) 
(5 Id) 



■K 4 - 



a: 



- (a) 




= 




!/:•! 

T/a=1 

T/co=3 
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300 
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FIG. 4: (Color online) (a) Plot of the relative fluctuation (Az)^_/z^ 
given in Eq. J52b vs the scaled evolution time f = tot for different 
temperatures: T/cj = (dotted line), TJ(o = 1 (dashed line), T/cj = 
3 (solid line), (b) Amplificatory plot of figure (a) from £ = to 20. 
Other parameters are set as y/co = 0.01, and g'/to = 0.01 15. 



Then the steady-state fluctuation is 

(2n, A + l)[y 2 + 4tu k (to k -2g k )] 



(Az(~)) 2 



(52) 



4 r 2 + 4^-i6g 2 

In Fig. [5] we plot the steady-state relative fluctuation 
(Az(oo)) 2 1 /z 2 as a function of the temperature T. It can be 
seen from Fig. [5] that the steady-state relative fluctuation in- 
creases with the increase of the temperature T. We can see a 
transition from squeezing to nonsqueezing when the temper- 
ature across a critical temperature T c , which can be obtained 
from Eq. (l52l for the case of k — — 1, 



T 

CO 



1 



T-r + % + 



3) 



(53) 



According to the parameters, we calculate the critical temper- 
ature in Fig.[5]is T c /a> = 0.222. 



V. DISCUSSIONS AND CONCLUSIONS 

We note that the two types of subsystems in our scheme, the 
superconducting charge qubit and the cantilever, have been 
well prepared in current experiments. Hence, it is possible 
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T/m 

FIG. 5: (Color online) Plot of the steady-state relative fluctuation 
(Az(oo)) 2 ^ 1 /^ against the scaled temperature T/a>. Other parameters 
are set as those in Fig. [4] Figure shows a transition from squeezing 
to nonsqueezing with the increase of the temperature of the environ- 
ment. 



to experimentally realize the scheme proposed in this paper 
within the reach of present-day techniques. Here we give 
a possible estimation of coupled-system parameters based 
on these published experimental parameters of the charge 
qubit and the cantilevers. The most important parameters in 
our model are the two coupling strengths g and g' given in 
Eqs. ( TTZb and ( |34j >, respectively. According to the current ex- 
perimental conditions dlz2]> we take the following param- 
eters. As an example, we choose a cantilever with a fun- 
damental frequency a> x 2n x 2 MHz, y = 2n x 20 kHz, 
and zo ~ 5 x 10~ 13 m. A magnetic tip produces a mag- 
netic gradient of \dB z /dz\ ~ 10~ 7 T/m at an approximate dis- 
tance 50 nm above the superconducting loop. Then we obtian 
C = \dB z /dz\zo ~ 5 x 10~ 6 T. We choose the area of the super- 
conducting loop S ~ 10~ 12 m 2 , and Ej « 5 GHz. Then we get 
4> ~ 2 An X 10~ 3 , which is suitable for making the approxima- 
tions of expanding the sine and cosine functions up to the first 
and second orders, respectively. Accordingly, the coupling 
strengths g * -2n x 6 MHz and g' * 2n x 23 kHz. Notice that 
the figures in the above sections are plotted in terms of these 
parameters. 

In conclusion, we have designed a theoretical scheme to re- 
alize tailorable couplings between a cantilever and a supercon- 
ducting charge qubit. By choosing proper parameters, both 
linear and nonlinear couplings can be achieved. We have also 
shown how to generate superposed coherent states and dy- 
namical squeezing in the cantilever based on the obtained cou- 
plings. We have investgated the influence of the environment 
on quantum states of the cantilever. It has been indicated that 
decoherence induced by the environment can drive the can- 
tilever from superposed coherent states into the steady coher- 
ent state. When the cantilever is initially in a coherent state, 
we have shown that there exists periodic position squeezing 
for the cantilever. Especially, it is found that under the action 
of the environment the cantilever can evolve from the vac- 
uum state to a steady state with position squeezing under a 



critical temperature of the environment. Therefore, the envi- 
ronment can induce the steady-state position squeezing of the 
cantilever. This reveals a new mechanism to create the steady- 
state squeezing and sheds new light on production of nonclas- 
sical effects of the cantilever. Finally, it should be emphasized 
that the experimental realization of the scheme proposed in 
the present paper deserves further investigation. 
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Appendix A: Derivation of quantum master equation (26l 

In this Appendix, we give a detailed derivation of the trans- 
form from quantum master equation ( l22b to equation d26*b . 
Starting from the quantum master equation 

p — i[p, a>a^a + g s (a^ + a)] 

y - t t t 

+— (n + l)(2apa' — a' ap — pa' a) 

7- iff 

+—n(2a pa — aa'p-paa ), (Al) 
we first make a rotation transform, 

p m = (0) p r (o) , r (0) = e - ie " ta , (A2) 

then we obtain p=R (0)p a) R i (0) and 

p = -i0a i aR(0)p m R 1 ' (0) + R (0) p m R* (0) 

+WR(0)p m a Jf aRH6). (A3) 

Substitution of Eq. ( IA3t into quantum master equation ( IAU 
leads to 

p (1) = i(a)-ff)p m a f a-i(co-b)a<ap m 

+ig sP (l) (aV' e + ae- w ) - ig s (aV e + ae- i6 )p m 

+| (n + l)(2ap (1) fl t - a f ap m -p^a^a) 

+ |n (2a t p (I) fl - flflV 1 ' - p (1 W) , (A4) 

where we have used the relations 

[R (60 , a 1 ' a] = [fl f (0) , a f a] = 0, (A5a) 

(#) a l R (0) = e i&,tfl a V iftjtfl = a V fl , (A5b) 
R 1 (6>) aR (0) = e l6a ' a ae-' l6a " a = ae^ 6 . (A5c) 

We choose a proper to ensure cj - 6 = 0. Under the initial 
condition 0(0) = 0, we get 

0(f) = cjt. (A6) 
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Then we obtain 

p (1) = ig s [pV,aU i0)t +ae- iM ] 

+ |(n + 1) (2ap (1) a t - a*ap m - p (1) a + a) 

+|n (2a 1 p (1) fl - flflV" - p (1 W) . (A7) 

Since the first term at the right-hand side of Eq. ( IA7t is a driv- 
ing term, in the following, we make a displacement transform 

p (2) = D + (a)p (I) D(a), D(a) = e aa ^ a " a ^ (A8) 

then we obtain p (1) = D (a) p (2) D ' (») and 

p (1) = b(a)p i2) DUfl) + D{a)p {2) DUa) 

+D(a)p {1) {t)bHa). (A9) 

Here we need to calculate the expressions for D(a) and 
(a). Making use of D (a) — e~T~e tra e~ a " a , we have 

at \ J at x ' 

= ~(aa + aa)D(a) + aD(a)e a " a a* e - a " a 
-a*D(a)a 

= — — (ad* — da*)D(a) + D(a)(aa^ — a* a), 

(A10) 

where we have used the formula e a ' a cfie~ a ' a = a' + a*. Simi- 
larly, we obtain 

£)'(a) = — — (aa* — aa*)D\a) + (a*a — da')D'(a). 

(All) 

Substitution of Eqs. dA9l . ( lAlOb . and (lAllt into Eq. (tA7t . 
we obtain 

p (2) = | (S + 1) ( 2fl p( 2 ) fl t - « ' flp (2) _ p <2) fl t fl ) 

+^(2aV 2 ^-oa t p (2) -p (2 W) 

+ (or* - fe,*-** + |a«) (ap (2) (r) - pfa) 

- (d + foe** + |a) (flV 2) - p (2) a f ) . (A12) 

Therefore, if we choose a proper a such that 

a + ig s e ih " + J -a = 0, (A13) 

then we obtain the quantum master equation d26l ). The solu- 
tion of Eq. dAl 31 ) can be obtained as 

«(0 - T^H*-?-***) (A14) 

under the initial condition a'(O) = 0. 



Appendix B: Derivation of Eq. ( l29t 

In this Appendix, we derive in detail the evolution of the 
cantilever governed by quantum master equation ( l22l i. For an 
initial state p(0), the state of the cantilever at time t can be 
obtained through the processes p(0) — > p (2) (0) — > p (2) (f) — > 
p(f)- The relationship between states p(0) and p <2) (0) is 

p®(0) = D t (a(O))/? t (0(O))p(O)/?(0(O))D(a(O)) 

= P(0). (Bl) 

For the initial state 

h«O)> = %(|0> + «*|-#), (B2) 

we have 

p (2) (0) = A/£(|0 + ^|-#)(<)B| + e-*<-0). (B3) 

The evolution process from state p {1) (0) to p (2) (?) is gov- 
erned by the quantum master equation 

p (2) (f) = J -(n+ l)(2ap (2) (f)fl t -fl t flp (2) (f)-p (2) (f)a 1 fl) 

+|n (2a 1 p (2) fl - ofl+p® (0 - p^^Oflfl 1 ") ■ (B4) 

In principle, the above master equation can be solved by the 
superoperator method 17011 . However, for simplicity, we only 
consider the zero temperature case in the following. At zero 
temperature, the master equation reduces to 

p (2) = |(2ap (2) (f)a t -a 1 ap (2) (f)-p <2) (f)a 1 fl)- (B5) 
Denoting the two superoperators 

Jp (2 \t) = ap i2 \t)a\ (B6a) 
Lp (2 \t) = -- (a'ap (2 \i) + p {2 \t)a'a) , (B6b) 

then the map from initial state p (2) (0) to p (2) (0 at time t is 
determined by 

p (2) (f) = exp(yfL) exp[(l - exp(- r f))7]p <2) (0). (B7) 

According to the initial state given in Eq. dB3l l, we have 

p (2) (0 = Nl fye$) (j8e-¥| + e^e h ' \pe~%) (-^~?| 

+e^e K \-f3e^)(/3e-i\ + \ -fie' S ) ( -pe~ ? |) , 

(B8) 

where we introduce the parameter A ( - = 2^8| 2 (e~ y ' - 1). During 
the derivation of the above equation dB8t , we have used the 
formulas 

exp[x/]|r)(s| = exp[jcrs*]|r)(s[, (B9a) 
exp(ytL)\r)(s\ = exp[-(|r| 2 + | 5 | 2 )(1 - e"*)/2] 

x \re~?}(se~% , (B9b) 
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for coherent states |r) and \ s). 

For obtaining the state in the Schrodinger picture at time t, 
we use the relation p (1) (f) = D(a(t))p {2 \t)D' ] ' (a(t)) to obtain 

D(a(t))p m (t)rt(a(t)) 
= N 2 {\ [ a(t)+/3e-^)(a(t)+/3e-^\ 

+ \a(f) ~pe-$)(a(i)~pe-%\ 

+c'¥ riAi |a(r) -pe-*)(a{i) +^e"?|) (BIO) 

where we introduce the parameter A, = 
2lm(a(t)f3* exp(-yf/2)). Here we have used the formula 



for coherent state \/3). And then we use the relation p(f) = 
R(8)p (1) (t)RH6) to obtain the state 



+e i> e A,-,A, p_ } </g+| + |^_ > f (B12) 

where we introduce the following two parameters 



p\ = (or(?)±j8e" ¥ )e" 



(B13) 



D(a)|^) = e 



/Im(ci' J S* ) I 



(Bll) 
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